The dynamical spin structure factor for the anisotropic spin-1/2 Heisenberg chain 
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The longitudinal spin structure factor for the XXZ-chain at small wave-vector q is obtained using 
Bethe Ansatz, field theory methods and the Density Matrix Renormalization Group. It consists of 
a peak with peculiar, non-Lorentzian shape and a high-frequency tail. We show that the width 
of the peak is proportional to q^ for finite magnetic field compared to q'^ for zero field. For the 
tail we derive an analytic formula without any adjustable parameters and demonstrate that the 
integrability of the model directly affects the lineshape. 
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One of the seminal models in the field of strong correla- 
tion effects is the antiferromagnetic spin-1/2 XXZ-chain 
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where J > is the coupling constant and h a magnetic 
field. The parameter A describes an exchange anisotropy 
and the model is critical for —1 < A < 1. Recently, 
much interest has focused on understanding its dynam- 
ics, in particular, the spin [1] and the heat conductivity 
[2], both at wave- vector q = 0. A related important ques- 
tion refers to dynamical correlation functions at small 
but nonzero q, in particular the dynamical spin struc- 
ture factors fi — x,y,z [3]. These quantities 
are in principle directly accessible by inelastic neutron 
scattering. Furthermore, they are important to resolve 
the question of ballistic versus diffusive transport raised 
by recent experiments [4] and would also be useful for 
studying Coulomb drag for two quantum wires [5]. 

In this letter we study the lineshape of the longitudi- 
nal structure factor S^^{q,Lij) at zero temperature in the 
limit of small q. Our main results can be summarized 
as follows: By calculating the form factors F{q,uj) = 
(O a) (here |0) is the ground state and \a) an ex- 
cited state) for finite chains based on a numerical evalu- 
ation of exact Bethe Ansatz (BA) expressions [6, 7] we 
establish that S^^{q,Lu) consists of a peak with peculiar, 
non-Lorentzian shape centered at lu ^ vq, where v is 
the spin-wave velocity, and a high-frequency tail. We 
find that |_F((7,a;)| is a rapidly decreasing function of the 
number of particles involved in the excitation. In par- 
ticular, we find for all A that the peak is completely 
dominated by two-particle (single particle-hole) and the 
tail by four-particle states (denoted by 2p and Ap states, 
respectively). Including up to eight-particle as well as 
bound states we verify using Density Matrix Renormal- 
ization Group (DMRG) that the sum rules are fulfilled 
with high accuracy corroborating our numerical results. 



By solving the BA equations for small A and infinite sys- 
tem size analytically we show that the width of the peak 
scales like q^ for h ^ 0. Furthermore, we calculate the 
high-frequency tail analytically based on a parameter- 
free effective bosonic Hamiltonian. We demonstrate that 
our analytical results for the linewidth and the tail are 
in excellent agreement with our numerical data. 

For a chain of length N the longitudinal dynamical 
structure factor is defined by 

^''('?''^) - nT. e-^'^'-'') / dte-*(5|(i)5|,(0)) 
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Here = J2j '''''' ^^'^ \^) is an eigenstate with 
energy Ea above the ground state energy. For a finite 
system, S^^ {q, w) at fixed g is a sum of (5-peaks at the 
energies of the eigenstates. In the thermodynamic limit 
N — > oo, the spectrum is continuous and S^^ (q,^) be- 
comes a smooth function of lo and q. By linearizing the 
dispersion around the Fermi points and representing the 
fermionic operators in terms of bosonic ones the Hamil- 
tonian (1) at low energies becomes equivalent to the Lut- 
tinger model [8]. For this free boson model S^^'{q,uj) can 
be easily calculated and is given by 



S^^ {q,u;)^K\q\5(u;^v\q\) 



(3) 



where K is the Luttinger parameter. This result is a 
consequence of Lorentz invariance: a single boson with 
momentum \q\ always carries energy lo = v\q\, leading to 
a ^-function peak at this level of approximation. 

We expect the simple result (3) to be modified in vari- 
ous ways. First of all, the peak at a; ~ should acquire 
a finite width 7^. The latter can be easily calculated for 
the XX point, A = 0, where the model is equivalent to 
non-interacting spinless fermions. In this case the only 
states that couple to the ground state via are those 
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containing a single particle-hole excitation {2p states). 
As a result, the exact S'^^{q,uj) is finite only within the 
boundaries of the 2p continuum. For h ^ 0, one finds 
Iq ~ q'^/m for small q, where m = (Jcosfci?)"^ is the 
effective mass at the Fermi momentum kp. For h = 0, 
m^^ and the width becomes instead 7, « Jq'^/8. 
In both cases the non-zero linewidth is associated with 
the band curvature at the Fermi level and sets a finite 
lifetime for the bosons in the Luttinger model. Different 
attempts to calculate 7, for A ^ have focused on per- 
turbation theory in the band curvature terms [9] or in 
the interaction A [10, 11] and contradictory results were 
found. All these approaches have to face the breakdown 
of perturbation theory near lo vq. 

Since perturbative approaches show divergences on 
shell, our discussion about the broadening of the peak 
is based on the BA solution. The BA allows us to cal- 
culate the energy of an eigenstate exactly from a system 
of coupled non-linear equations [12]. For A = these 
equations decouple, the structure factor is determined 
by 2p states only and one recovers the free fermion solu- 
tion. For I A| <C 1 the most important excitations are still 
of the 2p type and one can obtain the energies of these 
eigenstates analytically in the thermodynamic limit by 
expanding the BA equations in lowest order in A. For 
h {i.e., finite magnetization s = {Sj)) this leads to 

/ 2A \ Q 

7„ = 4 J 1 -I sin /ci? cos kp sin^ - « . (4) 

" V J 2 m* ^ ^ 

for the 2p type excitations. We therefore conclude that 
the interaction does not change the scaling of 7, com- 
pared to the free fermion case but rather induces a renor- 
malization of the mass given by m ^ m* = m/{l + 
2AamkF /t!')- We have verified our analytical small A 
result by calculating the form factors numerically [7] . For 
all A, we find that excitations involving more than two 
particles have negligible spectral weight in the peak re- 
gion. In Fig. 1 we therefore show only the form factors 
for the 2p states and a typical set of parameters. The 
form factors, for different chain lengths, N, collapse onto 
a single curve determining the lineshape of S'^^{q,uj) ex- 
cept for a high-frequency tail discussed later. The form 
factors are enhanced near the lower threshold ijJl{q) a-nd 
suppressed near the upper threshold u)u{q) in contrast 
to the almost flat distribution for A = 0. The lineshape 
agrees qualitatively with the recent result in [13] pre- 
dicting a power-law singularity at ujL{q) with a q- and 
A-dependent exponent. The inset of Fig. 1 provides a 
numerical confirmation that 7^ ^ q^, with the correct 
pre-factor as predicted in (4) for kp = 27r/5. For zero 
field, the bounds of the 2p continuum are known analyt- 
ically [14] and lead to a scaling 7, ~ q'^ for — 1 < A < 1. 
Furthermore, for /i = and A = 1, an exact result for the 
2p contributions to the structure factor has been derived 
[15]. 
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FIG. 1: (Color online) Form factors squared for the 2p excita- 
tions and diflterent iV at A = 0.25, s = -0.1 and q = 27r/25. 
The inset shows the scaling of the width 7,. The points are 
obtained by an extrapolation A'' cxd of the numerical data. 
The solid line is the prediction (4), 7g — 0.356 g^. 



Calculating a small number of form factors for finite 
chains poses two important questions: 1) For finite chains 
S^^{q,uj) at small q is dominated by 2p excitations. Is 
this still true in the thermodynamic limit? 2) How much 
of the spectral weight does the relatively small number 
of form factors calculated account for? We can shed 
some light on these questions by considering the sum rule 
I{q) = {2tt)-^ ^ dujS''{q,uj) = iV-i(S'^5'ig) where the 
static correlation function can be obtained by DMRG. 
As example, we consider again A — 0.25, s = —0.1, 
q = 27r/25 with TV = 200. For this case we have calcu- 
lated 2, 220, 000 form factors including up to 8p excited 
states as well as bound states. Note, however, that this 
is still small compared to a total number of states of 2^*^°. 
In the DMRG up to 2400 states were kept and the ordi- 
nary two-site method was utilized but with corrections to 
the density matrix to ensure good convergence with pe- 
riodic boundary conditions [16]. The typical truncation 
error was then ~ 10^^° and within the accuracy of the 
DMRG calculation (3 parts in 10^) the 2, 220, 000 form 
factors account for 100% of I{q). 99.97% of the spectral 
weight is concentrated in hiq), the contribution caused 
by the qN/2TT = 8 single particle-hole type excitations 
at Lu ^ vq. With increasing N we observe an extremely 
slow decrease in hiq)', however, even for a system of 2400 
sites, /2(g) is only reduced by 0.13% compared to the 
N — 200 case. While this large N behavior definitely 
requires further investigation it may not be very relevant 
to experiments, where effective chain lengths are limited 
by defects. 

Another feature missed in (3) is the small spectral 
weight extending up to high frequencies a; ~ J. This 
is relevant in the context of drag resistance in quan- 
tum wires because of the equivalence of S^^{q, uj) and the 
density-density correlation function for spinless fermions 
[5]. To calculate the high-frequency tail we start from 
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the Luttinger model 



The tail of S^^{q, uj) for /i ^ is then given by 



(5) 



Here, ^(x) is a bosonic field and n(a;) its conjugated 
momentum satisfying [(j){x),U{y)] — iS{x — y). The 
slowly varying part of the spin operator is expressed as 
~ ^JK/'K^x4>■ Note that both v and K depend on A 
and h. In the language of the Luttinger model, the spec- 
tral weight at high frequencies is made possible by boson- 
boson interactions. Therefore, we add to the model (5) 
the following terms 



5n{x) = j?_ 
+ n+ 



+ A cos (A^fKKcj) + Akpx^ , 



(6) 



where (j)R.L are the right- and left-moving components of 
the bosonic field with cj) = {(pL — 0/j)/\/2- They obey 
the commutation relations [dx4>L,R{x),dx<j)L,R{y)] = 
^idyS{x — y). These are the leading irrelevant opera- 
tors stemming from band curvature and the interaction 
part. The amplitudes r]±, (±, Cs and A are functions 
of A and h. For h ^ the A-term (Umklapp term) is 
oscillating and can therefore be omitted at low energies. 
Besides, the (^-terms have a higher scaling dimension than 
the ?7± -terms, so the latter yield the leading corrections. 
For h = 0, on the other hand, particle-holc symmetry 
dictates that r]± = and we must consider the ^-terms 
as well as the Umklapp term. For 7^ <^ oj — v\q\ <C J 
it is safe to use finite order perturbation theory in these 
irrelevant terms. 

In the finite field case the tail is due to the r]+- 
interaction. This allows for intermediate states with 
one right- and one left-moving boson, which together 
can carry small momentum but high energy uj » v\q\- 
It is convenient to write the structure factor defined 
in (2) as S'"'{q,uj) = — 2 Imx'"**(9, ijj) where x^'^* = 
~{K/'!T){dx(j)dx4') is the retarded spin-spin correlation 
function. The correction at lowest order in 7/+ to the 
free boson result then reads 



5x (9, i^n) = D^R {q, iw„) -I- £>^"^ {q, i(j„) ILrl {q, i(J„) 



(0) 



Url {q, ioJn) = - 



(7) 



xD^°^ (x,r)Z)r {x,r). 



(0) 



where D^^^^{q,iuJn) = {dx<t>L,Rdx(pL,R) = ±\q\/{iuJnT 
v\q\) are the free boson propagators for the right- and left- 
movers, respectively, and Url {q,ioJn) is the self-energy. 



5S;i{q,u;) 



Kr^U^ 6{u;-v\q\) 
vn ijp' — v^q^ 



(8) 



For h ~Q & connection between the intcgrability of the 
XXZ-model and the parameters in the corresponding 
low-energy effective theory exists [17]. The intcgrability 
is related to an infinite set of conserved quantities where 
the first nontrivial one is the energy current defined by 
Jb = / dx3E{x) with dx3E{x) = i\H{x),l dyH{y)] [18]. 
For the Hamiltonian (6) we find 



JE 



[{dxHf - {dx'l>Rf] - 4C- [{dx<pLf - {dxCi>Rf] 



+ 2C3 [dxMdx<pRf - dxMdxcj^L?] + 



(9) 



where the neglected terms contain more than four deriva- 
tives. Now conservation of the energy current, [Je, H] = 
0, implies = [22]. The spectral weight at high fre- 
quencies is therefore given by the ^+ and A-terms only. 

The perturbation theory for the C-i--term is analogous 
to the one for the ry^-tcrm. Now the incoming left (right) 
boson can decay into one left (right) and two right (left) 
bosons. This contribution is then given by 

5S^^l {q,u;) = ^^Q^i^^' - v\^)6{uj - v\q\) . (10) 

For the Umklapp term, we calculate the correlations fol- 
lowing [19] and find 

5Sr {q,u;) = Aq\w^ - v^q^-^iuj - v\q\) , (11) 

where A = 8Tr^\'^K^{2vf-^^ /V^iAK). We remark that, 
in a more general non-intcgrablc model, the (^3 term in 
Eq. (6) leads to an additional contribution to the tail 
which decreases with energy and becomes large near 
LU ^ vq. The increasing tail found in the integrablc case 
implies a non-monotonic behavior of S^^{q,Lu). Eqs. (8), 
(10) and (11) are valid in the thermodynamic limit. 
Wc can extend these results for finite systems and ex- 
press them in terms of the form factors appearing in 
(2). For a given momentum q = 2iTn/N, the form fac- 
tors generated by integer dimension operators as in (8) 
and (10) will then be situated at the discrete energies 

= 2-Kvl/N with I = n + 2,n + A, - ■ ■ . The form factors 
belonging to (11), on the other hand, will have energies 

= 2'Kv{l + AK) /N with Z = n, n -h 2, • • • . 

To compare our field theory results for the tail with BA 
data for the form factors we have to determine the a pri- 
ori unknown parameters in the effective Hamiltonian (6). 
In general, they can only be obtained in terms of a small- 
A expansion. To lowest order in A, Eq. (8) reduces to the 
weakly interacting result in [5, 11]. We also checked that 
in this limit Cs = for the XXZ-mode\ but Cs becomes 
finite if we introduce a next-nearest neighbor interaction 
that breaks intcgrability. For an integrablc model the 
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coupling constants can be determined by comparing ther- 
modynamic quantities accessible by BA and field theory. 
Lukyanov [20] used this procedure to find a closed form 
for and A in the case h = 0. Similarly, the parame- 
ters r]± can be related to the change in v and K when 
varying h and we find Jr]-{ho) = yj ^-k j Kv^ {a + 6/2)/6 
and Jr]^{hQ) — yj2T^ j Kv^hj \ where a — v^^ dv / dh\h^hg 
and b = K^^dK/dh\h^ii„. A numerical solution of the 
BA integral equations for u, K for infinite system size 
then allows us to fix r]± accurately for all anisotropics 
and fields so that the formulas for the tail do not contain 
any free parameters. The comparison with the form fac- 
tors computed by BA for finite and zero field is shown 
in Figs. 2 and 3, respectively. We note that the energies 
of the eigenstates are actually nondegenerate and spread 
around the energy levels predicted by field theory (see 
inset of Fig. 2). 
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FIG. 2; (Color online) Form factors (dots) obtained by BA 
compared to formula (8) (line) for A = 0.75, s = —0.1, TV — 
600 and q = 27r/50. The form factors for the exact eigenstates 
at a; « 2ttvI/N are added and represented as a single point. 
The inset shows each form factor separately. The number of 
states at each level agrees with a simple counting based on 
multiple particle-hole excitations created around the Fermi 
points. 

In summary, we have presented results for the line- 
shape of S^^{q,uj) for small q based on a numerical eval- 
uation of form factors for finite chains. We established a 
linewidth 7^ ^ for /i 7^ by solving the BA equations 
analytically for small A. In addition, wc showed that the 
spectral weight for frequencies jq <Ci^ — w|<z| ^ J is well 
described by the effective bosonic Hamiltonian. We pre- 
sented evidence that the hneshape of S^^{q,Lo) depends 
on the integrability of the model. This becomes manifest 
in the field theory approach by a fine tuning of coupling 
constants and the absence of certain irrelevant operators. 
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FIG. 3: (Color online) Sum of form factors at « 2nv{l + 
4:K)/N (dots) and at ci; ~ 2nvl/N (squares) obtained by BA 
for A = 0.25, s = 0, iV = 600 and q = 27r/50. The solid line 
corresponds to (11), the dotted line to (11) with finite size 
corrections included [21], and the dashed line to (10). 
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